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ABSTRACT 

Context. Recent radio observations have unveiled the existence of a number of radio pulsars with spin-down derived magnetic 
fields in the magnetar range. However, their observational properties appears to be more similar to classical radio pulsars than 
to magnetars. 

Aims. To shed light on this puzzle we first have to determine whether the spin-down derived magnetic field values for these 
radio pulsars are indeed representative of the actual neutron star magnetic field or if they are polluted, e.g. by the effects of a 
torque from a fallback disk. 

Methods. To investigate this possibility, we have performed deep IR (J, H, Ks bands) observations of one of these high magnetic 
field radio pulsars (PSR J1119-6127) with the ESO VLT to search for IR emission which can be associated with a disk. 
Results. No IR emission is detected from the pulsar position down to J ~24, H ~23 and Ks ~22. 

Conclusions. By comparing our flux upper limits with the predictions of fallback disk models, we have found that we can only 
exclude the presence of a disk with accretion rate M ^ 3 x 10^^ g s~^. This lower limit cannot rule out the presence of a 
substantial disk torque on the pulsar, which would then lead to overestimate the value of the magnetic field inferred from P and 
P. We have also compared the upper limit on the IR luminosity of PSR J1119-6127 with the IR luminosities of rotation-powered 
pulsars and magnetars. We found that, while magnetars are intrinsically more efficient IR emitters than rotation-powered pulsars, 
possibly because of their higher magnetic field, the relatively low IR emission efficiency of PSR J1119-6127 suggests that it is 
more similar to the latters than to the former. 
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seconds. Furthermore, the X-ray luminosity of both SGRs 
and AXPs largely exceeds their rotational energy {Lx ~ 
100 X E), while the rotational energy of radio pulsars can 
easily account for their X-ray emission {Lx ^ 0.001 x E; 
Becker & Triimper 1997; Possenti et al. 2002). The proper- 
ties of AXPs and SGRs are well explained by the magnetar 
model which interpretes these objects as isolated neutron 
stars with magnetic fields B ~ 10 ^"^"^^ G (hence dubbed 
magnetars), consistent with their observed spin-down with 
the usual vacuum dipole losses In the magnetar model, 
the X-ray luminosity is thought to be powered by the mag- 



^ Through the paper magnetic fields are computed (in G) 
using the relation B ^ 3.2 x 10^^ (P Pf^, where P is the NS 
period (s), P is its spin-down rate (s s~^). 



1. Introduction 

High-energy observations performed in the last two 
decades have unveiled the existence of a few unusual 
classes of neutron stars (NSs; see Popov 2007 for a re- 
view). The Anomalous X-Ray Pulsars (AXPs) and the 
Soft Gamma-Ray Repeaters (SGRs) are among them the 
most peculiar objects (see Woods & Thompson 2006 for a 
recent review) . At variance with the majority of the NSs 
known so far, i.e. the radio pulsars, they are typically ra- 
dio quiet but show X-ray pulsations at periods of a few 

Send offprint requests to: R. P. Mignani 
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netic field decay, while radio pulsations were believed to 
be suppressed by processes such as the photon splitting, 
which inhibit pair-production cascades in magnetic fields 
greater than the "quantum critical field" Be = 4.4x lO^'^ G 
(Baring & Harding 1998). 

This dichotomy between the two different pulsar 
classes - radio pulsars with B < Be on one side, and mag- 
netars with B > Be on the other - was shaken by the 
discovery of radio pulsars with magnetic fields above Be 
(Camilo et al. 2000). Despite having such high magnetic 
fields, although lower than those of the magnetars, these 
high-magnetic field radio pulsars (HBRPs) do not be- 
have according to any of the known magnetars templates. 
First of all, they are radio pulsars, while pulsed radio 
emission has been discovered so far only in the transient 
magnetar XTE J1810-197 (Camilo et al. 2006). Second, 
only two HBRPs, PSR J1119~6127 (Gonzalez & Safi-Harb 
2003) and PSR J 1718-3718 (Kaspi & McLaughhn 2005), 
have been detected in X-rays so far, with luminosities 
Lx ~ 10'^^"'^^ erg s~^ almost two orders of magnitude 
lower than those of the magnetars and smaller than their 
E. Finally, HBRPs do not show bursting emission, either 
in X-rays or in 7-rays, while AXPs and SGRs instead do. 
These differences might be explained assuming e.g., that 
HBRPs are dormant transients, that their lower X-ray lu- 
minosities are a consequence of their lower magnetic fields, 
or simply assuming that different evolutionary paths or 
stages account for the different phenomenologies. 

Of course, one alternative possibility is that the spin- 
derived magnetic field values of the HBRPs are unreliable 
because e.g., they are overestimated by the extra torque 
produced by a fossil disk formed out of residual matter 
from the supernova explosion. Fossil disks around isolated 
NSs have been invoked over the years to explain a large 
variey of phenomena (e.g. Michler & Dressier 1981; Lin 
et al. 1991; Phinney & Hansen 1993; Podsiadwolski 1993; 
Chatterjee et al. 2000; Alpar 2001; Menou et al. 2001; 
Blackman & Perna 2004; Cordes & Shannon 2006), and 
at least in the case of the AXP 4U 0142-1-61, recent Spitzer 
observations possibly revealed the presence of one of these 
disks (Wang et al. 2006). Thus, if HBRPs do have fos- 
sil disks, they should be detectable through deep, high- 
resolution IR observations. Since the IR luminosity of a 
hypothetical disk is expected to be larger for X-ray bright 
pulsars due to the flux contribution from the reprocess- 
ing of the X-ray radiation (Perna et al. 2000; Perna & 
Hernquist 2000) , the primary candidates are obviously the 
HBRPs detected in X-rays. 

In this work we report on the results of our recent 
deep IR observations of PSR J1119-6127. The pulsar was 
discovered in the Parkes multi-beam survey (Camilo et 
al. 2000) with period P = 407 ms and period derivative 
P ^ 4.022 X 10^^^ s s which give a characteristic age of 
~ 1600 years, a rotational energy loss E ~ 2.3 x 10'^^ erg 
s-S and a magnetic field B - 4.1 x lO" G. PSR J1119- 
6127 is also one of the very few pulsars with a measure 
of the braking index of 2.9 ± 0.1 (Camilo et al. 2000). X- 
ray emission was first detected with Chandra (Gonzalez 



& Safi-Harb 2003) which also revealed a compact pulsar 
wind nebula, while X-ray pulsations were discovered with 
XMM-Newton (Gonzalez et al. 2005). 
The structure of the paper is as follows: IR observations 
and results are described in §2, while comparisons with 
disk models and with IR observations of other isolated 
NSs are discussed in §3 and §4, respectively. 

2. IR Observations 

2.1. Observations Description 

IR observations of PSR Jl 119-6127 have been performed 
in Service Mode on January 25th, February 23rd, 24th and 
28th 2006 with NAos COnica (NACO), an adaptive op- 
tics (AO) imager and spectrometer mounted at the fourth 
Unit Telescope (UT4) of the VLT. In order to provide the 
best combination between angular resolution and sensitiv- 
ity, NACO has been operated with the S27 camera with a 
corresponding field of view of 28" x 28" and a pixel scale 
of 0'.'027. As a reference for the AO correction we have 
used the GSC-2 star S111230317098 {V = 13.7), located 
29'.'5 away from our target. Unfortunately, no suitable ref- 
erence star was available within the small NACO S27 field 
of view, which makes our AO correction not optimal and 
more sensitive on small scale fiuctuations of the atmo- 
spheric conditions. The Visual (VIS) dichroic element 
and wavefront sensor (4500 — 10000 A) have been used. 
Observations have been performed in the ESO Johnson 
J(A = 12650 1;AA = 2500 1), H{X = 16600 1; AA = 
3300 A) and K^{X = 21800 1; AA = 3500 1) filters. 

To allow for subtraction of the variable IR sky back- 
ground, each observation has been split in two sequences 
of short randomly dithered exposures with Detector 
Integration Times (DIT) of 20 and 55 s, and NDIT repeti- 
tions along each point of the dithering pattern (see Table 
[ij . This yields a total net integration time of about 2500 
s per band, per exposure sequence. For each exposure, 
the instrument readout mode has been selected accord- 
ing to the used DIT in order to minimize the read out 
noise. Owing to the expected faintness of the target, the 
DIT/NDIT combination has been modified after the Kg 
band observation of the first night to allow for a better 
signal-to-noise in the single exposures and to allow for a 
better hot pixels rejection. 

For all our observations, the seeing conditions were on 
average below 0"8 and the airmass was better than 1.3, al- 
lowing for a better yield of the NACO adaptive optics. Sky 
conditions were photometric in both nights. Night (twi- 
light flat fields) and day time calibration frames (darks, 
lamp flat fields) have been taken daily as part of the 
NACO calibration plan. Standard stars from the Persson 
et al. (1998) fields have been observed at the beginning of 
all nights for photometric calibration. As we expect the 
photometry errors to be dominated by the target's counts 
statistic rather than by the accuracy of the photometric 
calibration, we have not acquired photometric standard 
star fields prior to each exposure sequence. 



R. P. Mignani et al.: VLT observations of PSR J1119-6127 



3 



Table 1. Summary of the NACO J, H, X-band observations of the PSR J1119-6127 field with the number of exposure 
sequences, the total number of exposures per filter, the DIT and NDIT, the average seeing and airmass. 



yyyy-mm.dd 


Filter 


N 


exp 


DIT (s) 


NDIT 


Seeing (") 


Airmass 


2006.01.25 


Ks 


1 


8 


20 


16 


0.76 


1.26 


2006.02.23 


Ks 


2 


30 


55 


3 


0.66 


1.30 


2006.02.24 


H 


2 


30 


55 


3 


0.61 


1.29 




J 


2 


30 


55 


3 


0.88 


1.27 


2006.02.28 


J 


1 


15 


55 


3 


0.59 


1.32 



2.2. Data Reduction and Analysis 

The data have been processed using the ESO NACO 
pipeline]^ and the science images reduced with the pro- 
duced master dark and flat field frames. For each band, 
and for each night, single reduced science exposures 
have been combined to produce cosmic-ray free and sky- 
subtracted images. The photometric calibration pipeline 
yielded average zero points of 23.03±0.02 and 23.08±0.03 
(Kg) for January 25th and February 23rd, respectively, 
24.08 ± 0.04 (J) and 23.94 ± 0.04 (H) for February 24th, 
and 24.1 ± 0.05 (J) for February 28th. The data have 
been reduced independently using procedures run under 
the eclipse package ^ yielding qualitatively similar data. 

As a reference for the position of PSR Jl 119-6127 we 
have used its radio coordinates a( J2000)=ll'' 19™ 14.30", 
5(J2000)= -61° 27' 49'.'5, which have an accuracy of 0'.'2 
(Camilo et al. 2000). The astrometry on the NACO image 
have been computed using as a reference 7 stars selected 
from the 2MASS catalogue. The pixel coordinates of these 
stars (all non saturated and evenly distributed in the field) 
have been measured by gaussian fitting their intensity pro- 
files using the specific function of the GAIA (Graphical 
Astronomy and Image Analysis) tooj^ while the fit to the 
a, 6 reference frame has been performed using the Starlink 
package ASTROA/Q The rms of the astrometric solution 
turned out to be « 0'.'09 per coordinate. After accounting 
for the 0'.'2 average astrometric accuracy of ^MASSj^ the 
overall uncertainty to be attached to the position of our 
target is finally 0'.'3. 

2.3. Results 

Fig. [l] shows the Kg and H band images of the PSR 
Jl 119-6127 field with the computed pulsar radio position 
overlaid. No potential counterpart is detected at the ex- 
pected position, with the closest object being detected 
~ 1(7 away from the edge of the error circle. The same 
is true also for the J band image. We thus conclude that 
both the pulsar and its putative disk are undetected in 
each of the three observing bands down to estimated lim- 
iting magnitudes of J ^--^24, H ~23 and Ks ^ 22. At the 



http: / /www. eso.org/observing/dfo/quality/NACO/ 
http: / /www.eso.org/projects/aot/eclipse/ 
star-www.dur.ac.uk/ pdraper/gaia/gaia.html 
http: / /star-www. rl.ac.uk/Software/software. htm 
http: / /spider. ipac.caltech.edu/staff/hlm/2mass/ 



same time, no diffuse emission is recognized which can 
be possibly associated with the X-ray pulsar wind nebula 
detected by Chandra (Gonzalez & Safi-Harb 2003). 

3. Discussion 

3.1. Comparison with disk models 

We have used the derived IR fiux upper limits to con- 
strain the range of parameters that a hypothetical fos- 
sil disk around the pulsar could have. If a disk were in- 
deed present and interacting with the pulsar magneto- 
sphere, then, as mentioned in §1 and detailed below, the 
B field inferred from P and P could be largely overesti- 
mated. The torque exterted by a disk on the star mag- 
netosphere can be written as (e.g. Menou et al. 2001) 
■/disk — ~ —2MR?^^l, where M is the disk accre- 
tion rate, i?in is the disk inner radius, and f2 — 2-k / P 
is the angular frequency of the pulsar. The fact that PSR 
J1119-6127 is detected in radio implies that i?in cannot 
be smaller than the light cylinder radius R\c = c/fi (e.g. 
lUarionov & Sunyaev 1975). On the other hand, if the 
inner radius of the disk were outside of the light cylin- 
der, where the magnetic field lines are open, no efficient 
torque could operate. Therefore, in the following analysis 
we consider only the case i?in — R\c, which yields a torque 
■/disk — —2Mc^/fl or, equivalently, an energy loss (in mod- 
ulus) Edisk ~ 2Afc2 = 2 10^7 A//(10i6gs-i) erg s^^. 
Under these conditions, the total energy loss of the pulsar, 
accounting for both the dipole and the disk torque compo- 
nents, is given by £: = £^dip+£^disk ^ B^Cl^R'^ /6c^+2Mc^, 
where R is the radius of the star. Clearly, if -Edisk ^ £'dipi 
the value of B that is inferred from P and P, i.e. assuming 
E = E^ip, could be largely overestimated. In the case of 
PSR J1119-6127, a fallback disk with accretion rate of the 
order of 10^^ g could account for the entire energy loss 
of the pulsar {E - 2.3 x 10^^ erg s'^), even without the 
contribution of dipole losses, which would be the case for 
a very low magnetic field. We thus take M ~ 10"'^^ g s~^ 
to be the accretion rate corresponding to the maximum 
torque that could be produced by a hypothetical disk. 

We have simulated the disk IR spectrum using the disk 
model developed by Perna et al. (2000), which takes into 
account the contribution to the disk IR emission due to 
both viscous dissipation and reprocessing of the X-ray ra- 
diation from the pulsar. The spectra have been renormal- 
ized for the distance d to PSR J1119-6127. Camilo et al. 
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Fig. 1. 6" X 6" sections of the VLT/NACO and H band images of the PSR J1119-6127 field. North to the top, 
East to the left. The effect of the worse AO correction (see §2.1) is recognized from the asymmetric PSF of the stars in 
the field. The circle (O'.'S radius) corresponds to the pulsar radio position uncertainty after accounting for the accuracy 
of our astrometric solution (see §2.2). 




(2000) reported d = 2.4 - 8 kpc, while Gonzalez & Safi- 
Harb (2003), based on the measured extinction per unit 
distance in the pulsar direction, estimated d = 5.4 — 12.6 
kpc. Most likely, the pulsar is not further than 8 kpc, ac- 
cording to its location with respect the Carina spiral arm 
(Camilo et al. 2000). In the following we report our results 
as a function of DQ=d/ (6 kpc). Fig. 2 shows the modelled 
disk IR spectra computed for two different values of the 
disk accretion rate M compared to the observed IR flux 
upper limits. In particular, we show the predicted flux cor- 
responding to the maximum value of M that would make 
the disk emission compatible with the current limits, as 
well as the flux corresponding to the maximum value of 
M compatible with the spin down rate of the pulsar. 

We find that our limits only rule out disks with M ^ 
3x 10^^ Dg g s~^, i.e. well above the value of 10^^ g s~^ cor- 
responding to the case of maximum allowed disk torque. 
Therefore, we cannot exclude with the current observa- 
tions that the magnetic field derived from the pulsar spin- 
down, under the assumption of a purely vacuum dipole 
energy loss, be overestimated due to pollution by a disk 
torque. We note that, given the low X-ray luminosity of 
PSR J1119~6127, io.5-iokcV = S.Stg^g x 10^2 x I?| erg s^^ 
(Gonzalez & Safi-Harb 2003), the contribution to the disk 
IR emission due to the X-rays reprocessing is so low that 
it becomes comparable with that due to viscous dissipa- 
tion in the disk only for M ^ 5 x 10^^ D§ g s~^. For this 
reason, the derived IR flux upper limits for PSR J1119- 
6127 are less stringent in ruling out a fallback disk at the 
light cylinder with respect to similarly deep upper limits 
obtained for the AXPs which, instead, have a much higher 
X-ray luminosity. Indeed, an X-ray luminosity higher by 
a factor 100 would raise the disk IR emission much closer 
to our present upper limits. 

We have seeked for other evidence which might indi- 
rectly unveil the existence of a fallback disk and the effect 



Bio-' 



K H J 



M„=0,3 



1 I \ 



M„=0,01 



i>(Hz) 

Fig. 2. Expected IR spectrum of a fallback disk of in- 
ner radius i?i„ = R\c and two values of the disk ac- 
cretion rate M: solid line: maximum value of M com- 
patible with the IR limits; dashed line: maximum value 
of M compatible with the pulsar's spin down rate. The 
IR flux upper limits have been corrected for the inter- 
stellar extinction applying the relations of Fitzpatrick 
(1999) for Av— 5, as derived from the X-ray absorption 
Nh = 911 X 10^^ cm-2 (Gonzalez & Safl-Harb 2003) and 
the relation Ay — Nh /1-79 x 10^^ atoms cm"^ mag 
(Predhel & Schmitt 1995). 



of its torque on the pulsar's spin down. In principle, a 
torque from a fallback disk should leave a signature in 
the pulsar timing by increasing the level of the timing 
noise. In the case of PSR J1119-6127, the level of the ra- 
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dio timing noise does not show any clearly anomalous ex- 
cess which might be associated with the effect of an acting 
disk torque, and it seems apparently consistent with the 
level expected for its high P (Arzoumanian et al. 1994). 
However, we note that the magnitude of the effect would 
depend on the actual value of the disk torque, which is 
obviously unknown, and it might be confused with the 
underlying timing noise. In any case, the timing analysis 
can in no way rule out that the pulsar's spin down might 
have been affected by a disk torque in the past. 

3.2. Comparison with otiner NSs 

For the estimated ranges of distance and Nh (see previous 
section), our flux upper limits yield for PSR J 1119-6127 
a JsT-band IR luminosity Lk < 6.6l" x 10^° x Dj erg 
s~^. We have compared this upper limit with the IR lu- 
minosities of different classes of isolated NSs. We caveat 
here that the nature of the IR emission may be differ- 
ent across the whole sample. For instance, in the case of 
rotation-powered pulsars, the IR emission is thought to 
be produced in the NS magnetosphere, as shown by their 
power-law spectra (e.g. Shibanov et al. 2003), while in the 
case of the magnet ars it might be produced by an X-ray ir- 
radiated, (though passive), fallback disk, by the magnetic 
field decay, or by curvature radiation in the magnetars' 
coronae (Beloborodov & Thomspon 2007). Nevertheless, 
comparing the IR properties of different classes of isolated 
NSs can still be useful to unveil similarities and diversi- 
ties which can be indeed ascribed to different emission 
processes and thus be used to disantangle, e.g. magneto- 
spheric and disk emitters. 

Table |2] summarizes the IR flux measurements for all 
the isolated NSs with an IR counterpart, i.e. rotation- 
powered pulsars and magnetars. In order to make a con- 
sistent comparison with rotation-powered pulsars, which 
are persistent emitters, for the magnetars we have selected 
only IR flux measurements taken when the X-ray source 
was as close as possible to quiescence. We have include 
in our compilation also the AXP IRXS J170849-400910, 
although its IR identification has not been confirmed yet 
(Safi-Harb & West 2005; Durant & van Kerkwijk 2006; 
Rea et al. 2007a), hence we did not consider it in the fol- 
lowing analysis. The proposed identification of IE 1841- 
045 (Wachter et al. 2004) has been discarded by high- 
resolution IR observations (Durant 2005). No IR emis- 
sion has been detected so far from the X-ray Dim Isolated 
NSs (XDINSs; Mignani et al. 2007; Lo Curto et al. 2007; 
Rea et al. 2007b) and from any compact central objects 
(CCOs) in SNRs (Wang, Kaplan & Chakrabarty 2007; 
Fesen, Pavlov & Sanwal 2006). 

For each object we have computed its IR luminosity ei- 
ther in the Ks or in the K band, as we estimate the error 
due to the passband correction to be below ^ 0.1 mag- 
nitudes, i.e. fully acceptable for the goals of our analysis. 
For Vela and Geminga, we have extrapolated their if-band 
magnitudes from the IR colors. Passband transformations 



between different if-band filters have been neglected. The 
flux conversion from the HST/ NICMOS passbands to the 
Johnson's ones has been applied using the NICMOS Units 
Conversion tooQ 

For the rotation-powered pulsars, distance values have 
been taken either from the available radio/optical par- 
allaxes or from the radio dispersion measure (Cordes & 
Lazio 2002)0 For the magnetars we have used either the 
distances of the parental stellar clusters or of the associ- 
ated supernova remnants, or the distances inferred from 
the Nh (see Table |2] and references therein). For the inter- 
stellar extinction correction we have applied the relations 
of Fitzpatrick (1999) using as a reference either the mea- 
sured Ay or the value derived from the Nj^ recomputed 
from our X-ray spectral fits and the relation of Predehl 
& Schmitt (1995). For the magnetars we have fitted an 
absorbed power-law plus a blackbody model (see Tiengo 
et al. 2005; Rea et al. 2004; Rea et al. 2005; Woods et 
al. 2004; Patel et al. 2003; Morii et al. 2003; Mereghetti 
et al. 2004 for further details on the single observations) 
over the spectral range 2-10 keV. All the Nh values have 
been computed assuming solar abundances from Anders 
& Grevesse (1989). Although the reference Ay have been 
obtained with different methods, this does not affect sig- 
nificantly our estimates of the IR luminosity, expecially 
in the i^T-band where the effects of the interstellar extinc- 
tion are lower. The overall IR luminosity errors take into 
account the measured photometric errors as well as all 
the uncertainties on the isolated NS distance and on the 
interstellar extinction correction, all reported in Table [2j 

3.3. Results 

In the left panel of Fig. [3]we have plotted the computed IR 
luminosities Lk for all the isolated NSs in Table [2] and the 
upper limit for PSR J1119-6127, as a function of the NSs 
rotational energy loss E. From this plot we clearly see that 
the rotation-powered pulsars and the magnetars cluster in 
quite distinct regions of the diagram. In particular, PSR 
J1119-6127 is definitely closer to the group of the rotation- 
powered pulsars than to the magnetars one, which would 
suggest its connection with the formers rather than with 
the latters. 

From a general point of view, it is clear that there is 
no substantial difference between the IR luminosity of the 
magnetars and that of the young (< 5, 000 years) rotation- 
powered pulsars (Crab and PSR B1509-58), which have 
all luminosities Lk ~ 10^° — lO'^^ erg to be compared 
with Lk ~ 10^^ — 10^^ erg s~^ of the older rotation- 
powered pulsars (Vela, PSR B0656+14 and Geminga). 
For the rotation-powered pulsars Fig. |3] shows that the 
IR luminosity correlates rather well with the rotational 
energy loss, with Lk oc e^-3±0-04^ This correlation is 
similar to the one found for the optical luminosity, i.e. 
Lopt oc £;i-6±o.2 ^ggg^ g Kramer 2004), which confirms 

^ http: / /www. stsci.edu/hst/nicmos /tools / 

® http://rsd-www.nrl.navy.mil/7213/Iazio/ne_model/ 
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Table 2. Summary of the IR fluxes measurements for all types of isolated NSs with an identified IR counterpart i.e. 
rotation-powered pulsars (rows 1-5) and magnetars (rows 6-11). The columns give the observed J, H, K, Kg magnitudes 
(an hyphen stands for non-detection, values in italics have been extrapolated), the distance, and the interstellar 
extinction Ay either derived from existing optical measurements (O) or from the Nh derived from the fits to the 
X-ray spectra (X) by using the relation of Predehl & Schmitt (1995). K-hanA flux values in italics have been derived 
from the extrapolation of the J and TJ-band fluxes. 



NS Name 


J 


H 


K 


Ks 


d(kpc) 


Av 


Ref. 


Crab 


14.8±0.05 


14.3±0.05 


13.8±0.05 




1.730±0.28 


1.62 (0) 


1,2,3 


PSR B1509-58 




20.6±0.20 




19.4±0.1 


4.181±0.60 


4.8 (0) 


4,2,5 


Vela 


22.7±0.10 


22.0±0.16 


21.3±04 




n 294+0 019 


0.20 (0) 


6,7,8 


PSR B0656+14* 


24.4±0.10 


23.2±0.08 


22.6±0.13 




U.^8»_0.027 


0.09±0.06(O) 


9,10,11 


Geminga* 


25.1±0.10 


24.3±0.10 


23.4±0.4 




n 11-7+0.059 
u.ro/ _o.034 


0.12±0.09(O) 


9,12,13 


4U 0142-^61** 






19.7±0.05 


20.1±0.08 


>5 


5.1 (X) 


14,15 


IE 1048.1-5937 


21.7±0.30 


20.8±0.30 




21.3±0.30 


3±1 


6.1 (X) 


16,17,15 


IRXS J170849-400910'' 


20.9±0.10 


18.6±0.10 




18.3±0.10 


5 


7.8 (X) 


18,15 


XTE J1810-197 




22.0±0.10 




20.8±0.10 


4±1+ 


5.1 (X) 


19,20,15 


IE 2259+586 








21.7±0.20 


3.0±0.5 


5.7±0.1 (O) 


21,22,23 


SGR 1806-20 


20.1±0.14 


15.1 ± 1.6 


29±2 (0) 


24,25,26 



(1) SoUerman (2003); (2) radio dispersion measure, Cordes & Lazio (2002); (3) Sollerman et al (2000); (4) Kaplan & Moon 
(2006); (5) Lortet et al. (1987); (6) Shibanov et al. (2003); (7) radio parallax, Dodson et al. (2003); (8) Mignani et al. (2003); 
(9) Koptsevich et al. (2001); (10) radio parallax, Brisken et al. (2003); (11) Pavlov et al. (1997); (12) optical parallax, Caraveo 
et al. (1996); (13) Kargaltsev et al. (2005); (14) HuUeman et al. (2004); (15) present work; (16) Wang & Chakrabarty (2002); 
(17) Gaensler et al. (2005); (18) Israel et al. (2003); (19) Israel et al. (2004); (20) Rea et al. (2004); (21) HuUeman et al. (2001); 
(22) Kothes et al. (2002); (23) Woods et al. (2004); (24) Israel et al. (2005); (25) Mc Clure-Griffiths & Gaensler (2005); (26) 
Eikenberry et al. (2004) 

* magnitudes refer to the HST/NICMOS filters HOW, 160W, 187W, which overlap the J, H and K passbands 
** K and Ks magnitudes have been taken at different epochs 

^ IR counterpart still to be confirmed (Safi-Harb & West 2005; Durant & van Kerkwijk 2006; Rea et al. 2007a) 
+ revised downward to 2.5 kpc (Gotthelf & Halpern 2005) 

that the IR emission of rotation-powered pulsars, like the 
optical one, is mostly magnetospheric. 

Instead, for the magnetars the scatter of the points 
does not allow to recognize a correlation between Lx and 
E. However, if the magnetars' IR emission were also pow- 
ered by their rotational energy they would be much more 
efficient IR emitters than the rotation-powered pulsars, 
with IR luminosities much closer to their E. In partic- 
ular, we note that, if the distance lower limit of 5 kpc is 
confirmed, the IR luminosity of 4U 0142-1-61 could be com- 
parable to its its E, making it the intrinsically more lumi- 
nous magnetar. This intrinsically larger IR output could 
be explained by the presence of an additional source of 
emission which, at least in the case of 4U 0142-1-61, might 
be identified with a possible fossil disk (Wang et al. 2006). 
The same might be true also for the other magnetar with 
the highest IR luminosity (SGR 1806-20), while for the 
others the presence of a surrounding fallback disk appears 
less compelling. Alternatively, it is possible that the IR 
emission of magnetars is powered, as it is in the X-rays, 
by the star magnetic field rather than by its rotation. We 
have plotted in the right panel of Fig. [3] the IR emission 
efficiency as a function of the dipole magnetic field B. 
Despite the scatter of the points, it is clear that the mag- 
netic field does imply a larger IR emission efficiency for the 
magnetars than for the rotation-powered pulsars. We thus 



speculate that, although the contribution of a disk cannot 
be a priori ruled out, the IR emission of the magnet axs 
is substantially driven by the magnetic fleld. In particu- 
lar, we note that with a magnetic field B ^ 4.1 x lO^'^ G, 
one might expect for PSR J1119-6127 a magnetar-like IR 
emission efficiency, while it is at least one order of magni- 
tude lower. This makes PSR J1119-6127, once again, more 
similar to the rotation-powered pulsars than to the mag- 
netars. This might suggest that the actual magnetic field 
of PSR Jl 119-6127 is lower than the measured one and 
that a torque from a disk might have indeed affected the 
pulsar's spin down. However, we note that, given the disk 
accretion rate compatible with the maximum torque and 
the low X-ray luminosity of the pulsar (see §3.1), the con- 
tribution of such a disk to the total IR flux would likely 
be low enough for the pulsar IR emission to be domi- 
nated by the magnetospheric component, as in the classi- 
cal, rotation-powered radio pulsars. 

4. Conclusions 

We have reported on deep IR observations performed with 
the ESO VLT to constrain the presence of a fallback disk 
around the high magnetic fleld radio pulsar PSR J 1119- 
6127. No IR emission has been detected at the pulsar's po- 
sition down to limiting magnitudes of J ^^24, H ^23 and 
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Fig. 3. left panel: measured if -band luminosities for all the isolated NSs listed in Table [2| as a function of the NSs 
rotational energy loss E. For J1119-6127 (labelled in the figures) we have plotted the IR luminosity upper limits 
corresponding to the most extreme distance/absorption values. For 4U 0142+61 we have plotted the IR luminos- 
ity lower limit corresponding to a lower limit on the source distance of 5 kpc (HuUemann et al. 2004). Red filled 
circles and blue filled squares indicate rotation-powered pulsars and magnetars, respectively. The red dashed line cor- 
responds to the linear fit for rotation-powered pulsars while the solid line shows the limit case Log{LK) = Log(E). 
Right panel: the derived IR eSiciencies, defined as the ratio between Lx and E, as a function of the dipole mag- 
netic field B inferred from the NS spin down. Timing parameters have been taken from the ATNF Pulsar Database 
(http:/ /www. atnf.csiro.au/research/pulsar/psrcat) for rotation-powered pulsars, and from Woods & Thompson (2006) 
for the magnetars. For both IE 1048-5937 and SGR 1806-20 we have taken as a reference the average P value. 



Ks ~22. These upper limits have been compared with the 
expected IR spectrum emitted from a fallback disk, which 
we have computed using the disk models of Perna et al. 
(2000). We have found that the current fiux upper limits 
only rule out a fallback disk with M ^ 3 x 10^^ g s^^. 
However, a disk with an accretion rate of ~ 10^^ g s~^ 
can still account for the rotation energy loss of the pulsar, 
hence we cannot yet confirm or exclude that the pulsar ex- 
periences an extra torque produced by a fallback disk, and 
that the value of the magnetic field inferred from P and P 
is thus overestimated. We have also compared the upper 
limit on the IR luminosity of PSR J1119-6127 with the 
measured IR luminosities of rotation-powered pulsars and 
magnetars. While magnetars are intrinsically more effi- 
cient IR emitters than rotation-powered pulsars, probably 
because of their higher magnetic field, we have found that 
the relatively low IR emission efficiency of PSR J1119- 
6127 makes it more similar to the latters than to the for- 
mers. Although not strictly compelling, this might be an 
indication of a magnetic field actually lower than the mea- 
sured one. 
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